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Abstract 

Let G = (V, E) be a planar n-vertex digraph. Consider the problem of computing max 
si-flow values in G from a fixed source s to all sinks t G V \ {s}. We show how to solve 
this problem in near-linear 0(n log 3 n) time. Previously, no better solution was known than 
running a single-source single-sink max flow algorithm n — 1 times, giving a total time bound 
of 0(n 2 logn) with the algorithm of Borradaile and Klein. 

An important implication is that all-pairs max si-flow values in G can be computed 
in near-quadratic time. This is close to optimal as the output size is 8(n 2 ). We give a 
quadratic lower bound on the number of distinct max flow values and an 51 (n 3 ) lower bound 
for the total size of all min cut-sets. This distinguishes the problem from the undirected 
case where the number of distinct max flow values is O(n). 

Previous to our result, no algorithm which could solve the all-pairs max flow values 
problem faster than the time of 0(n 2 ) max-flow computations for every planar digraph was 
known. 

This result is accompanied with a data structure that reports min cut-sets. For fixed s 
and all i, after 0(n 3 / 2 log 3 / 2 n) preprocessing time, it can report the set of arcs C crossing 
a min si-cut in time roughly proportional to the size of C. 
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1 Introduction 



Computing max flow in a graph is a classical algorithmic problem dating back to Ford and 
Fulkerson |13[ 114] . Given a graph G = (V,E) with arc capacities, a source s E V and a sink 
t G y\{s}, the problem is to send as much flow as possible from s to t without violating capacity 
constraints and flow conservation at vertices in V\{s,i}. Many polynomial-time algorithms for 
this problem are known. In this paper we focus on planar graphs. There exists an 0{n log n) time 
algorithm for planar digraphs |3j[T0] and an 0(n log log n) time algorithm for planar undirected 
graphs |21] , More recently, an 0(n log 3 n) time algorithm for a more general problem with a set 
of sources and a set of sinks was proposed [3]. 

It is natural to consider the problem of finding multiple max flows and min cuts, e.g., to 
study the connectivity of the graph. For undirected (not necessarily planar) graphs, min si-cuts 
for all source/sink pairs (s, t) form a laminar family and can be compactly represented by a tree 
structure known as a Gomory-Hu tree [16] . Finding a Gomory-Hu tree reduces to solving n — 1 
min si-cut problems |17] , For planar undirected graphs, this gives an 0(n 2 log log n) time bound 
using the algorithm in |21] , This time bound was significantly improved to 0(n log 5 n) in [6]. 
In particular, this means that all max si-flow/min si-cut values can be computed in near-linear 
time. 

We consider the same problem but for planar digraphs. This problem is more challenging 
since the Gomory-Hu property no longer holds, even for planar digraphs [2]. In fact, there is an 
example of (non-planar) digraphs with 0(ra 2 ) distinct max flow values [22]. We show that this 
bound holds for planar digraphs as well. Figure [1] shows a plane n-vertex digraph with 0(n 2 ) 
different cycles. There is a well-known duality between shortest cycles in the primal G and min 
si-cuts in the dual G* of a plane graph: let / and g be distinct faces in G and let /* and g* 
be the corresponding dual vertices in G* , respectively. Let C be a shortest clockwise cycle in G 
such that / is outside of C and g is inside. Then there is a min f*g*-cut in G* such that the 
arcs crossing this cut are exactly the (dual) arcs of C. It is now easy to see that the dual of the 
graph in Figure [1] has 0(n 2 ) distinct min cuts. Moreover, the total size of the min cut-sets in 
this example is 0(n 3 ). 

At first sight it might thus seem impossible to compute the capacities of all min cuts in 
o(n 3 ) time. Indeed, repeated applications of the single-source single-sink max flow algorithm of 
Borradaile and Klein [3] leads to a time bound of 0(n 3 logn). Arikati et al. |T] computed all of 
these capacities in 0(n 2 +7 3 log 7) time where 7 is the minimum number of hammocks, a special 
kind of outerplanar subgraphs, into which G can be partitioned. However, 7 might be O(n) and 
the running time of this algorithm is also 0(n 3 log n) in this case. Nevertheless, we are able to 
significantly improve this to 0(n 2 log 3 n) which is optimal up to logarithmic factors. In fact, we 
show something stronger: for fixed source s, max si-flow values for all n — 1 sinks i ^ s can be 
found in a total of 0(n log 3 n) time. 

Our algorithm can be changed to a data structure that can report the arcs crossing the min 
cuts, i.e., min cut-sets. Given a fixed source s it requires 0(n 3 / 2 log 3 / 2 n) preprocessing time. 
Afterwards, it can report a min cut-set C for s and any sink t ^ s in 0(|C| log log n) time. 
Hence, the total time to compute all min cut-sets C in the planar digraph is 0(n 5 / 2 log 3 / 2 n + 
SceC 1^1 log log n )- With a longer preprocessing of 0(n 5//3 log 2//3 n) time, we can decrease the 
report time to 0(|C|). 

Our result shows that in planar digraphs all 0(n 2 ) max flow values can be computed in time 
essentially equal to the time for finding a linear number of max flows. This joins the result of 
Arikati et al. [lj that showed an 0(n 2 ) algorithm for the same problem for bounded treewidth 
digraphs. For general n-vertex digraphs, Hao and Orlin |18] showed how to compute O(n) 
max flows in the time it takes for a single max flow computation and used this result to find a 
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Figure 1: Example of a planar digraph on 4n vertices for which min /*^-cuts in the dual are 
pairwise distinct, 2 < i,j < n. There are thus 0(n 2 ) min cuts and each has 0(n) arcs for a 
total of 0(n 3 ). In order to make the dual graph simple one can triangulate the above graph 
with infinite weight edges. 

(global) min cut. Based on this and our result for planar digraphs, we conjecture that also in 
general digraphs all 0(n 2 ) max si-flow values can be computed faster than the time required 
for computing 0(n 2 ) max si-flows separately. The result of Hao and Orlin does not resolve this 
conjecture as the 0(n) source/sink pairs considered by their algorithm are in a sense arbitrary; 
in particular, their result cannot be used to efficiently find max si-flow values for fixed s and all 
t / s. 

When all arc capacities are equal to a single unit, we get the arc connectivity problem. For 
this special case of our problem, Cheung et al. [9] showed an 0{m u ) all-pairs arc connectivity 
algorithm, for an m-arc digraph, where oj is the matrix multiplication exponent. This is faster 
than running the 0(min{m 1//2 , n 2//3 } • m) algorithm for si-arc connectivity of Even and Tarjan 
for every pair of vertices when m = 0(n L94 ), using the currently best known value of 
u> < 2.3727 |28j . Cheung et al. also showed an 0(d UJ ~ 2 n UJ ' 2+1 ) algorithm for the same problem 
for well-separable digraphs, which also include planar graphs, with maximum degree d. 

Cabello et al. [7] gave a simple algorithm for computing the shortest non-contractible cycle 
in a directed graph embedded in a sphere with b boundaries, which uses b minimum cut compu- 
tations with a fixed source in the dual planar graph. Our result improves the time bound of this 
algorithm from 0(bn log n) to 0(n log 3 n) when b = w(log 2 n). In a similar way, our algorithm 
can be used to find the shortest clockwise cycle in an embedded directed planar graph, by finding 
the min si-cut in the dual planar graph where s is the vertex dual to the inifinite face of the 
original graph. 

The organization of the paper is as follows. In Section [2j we make some definitions and 
simplifying assumptions. In Section[3l the overview of the algorithm for our problem is presented. 
Our algorithm is guided by a recursive decomposition of the input graph. It recursively computes 
max preflows to the outer boundaries of all subgraphs in the decomposition using previously 
computed max preflows to outer boundaries of the parent subgraphs. When the algorithm 
reaches a leaf of the recursive decomposition, a max preflow to a sink i contained in that leaf 
is found. The value of a max si-flow can then be identified as the amount of preflow into i. A 
straightforward implementation of this algorithm will not lead to a near-linear time bound. In 
Section [4j we show how to efficiently implement various steps of the algorithm. An important 
tool here is a modification of a flow fixing procedure from [3]. In Section [5] we show how the 
min si-cut-sets themselves can be reported. Finally, in Section [6] we give some conclusions and 
suggest future directions of research. 
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2 Preliminaries 



Let G = (V, E) be a simple planar digraph, where V is the vertex set and E C V x V is the 
set of arcs. Let n = \V\. Since G is planar, we have \E\ = 0(n). We assume that for every 
arc e = (u,v) 6 E the reversed arc rev(e) = (v,u) is also in E. This assumption can be easily 
satisfied for planar graphs, by using the same embedding for e and rev(e). By edge we mean a 
pair of two opposite arcs, e and rev(e), which share their embedding in the plane. 

We use the standard definition of the dual planar graph G* , that is, every dual vertex 
corresponds to a primal face, and two dual vertices are connected by an edge if and only if there 
is an edge between the two corresponding primal faces. We orient the dual arcs such that the 
dual arc of a primal arc e is oriented from the right side of e to its left side. If the primal is a 
flow network with capacities on arcs, the dual arcs have weights equal to these capacities. 

2.1 Flows 

For graph G the capacities of the arcs are given by a capacity function c : E — > M + . Let s £ V be 
a source and let t E V \ {s} be a sink. We define a flow assignment to be a function f : E —> M 
satisfying antisymmetry, i.e., for all e £ E we require /(e) = —f{rev{e)). A flow assignment is 
called an st-pseudoflow if for all arcs e £ E we have /(e) < c(e). The inflow of an st-pseudoflow 
/ for a vertex v G V is defined as: 

inflow/ (v) = ^ /(«,«), 

We say that st-pseudoflow / has excess in u when inflow/ (v) > and in this case we refer to 
inflow/ (v) as its excess (value). Similarly, we say that / has deficit in v when inflow/ (v) < and 
inflow f(v) is its deficit (value). An st-preflow is an st-pseudoflow such that no vertex in ^\{s} 
has deficit. An st-flow is an st-preflow such that no vertex in V \ {t} has excess. When no 
confusion arises, we shall omit s and t and simply talk about pseudoflows, preflows, and flows. 
A circulation is a pseudoflow such that for every vertex v £ V, inflow/ (y) = 0. A circulation 
does not have a source or a sink. 

The value of a preflow / is given as inflow/ (t). Finally, a max pseudoflow (preflow) is a 
pseudoflow (preflow) such that there is no residual path that starts at a source or at a vertex 
with excess and ends at a sink or at a vertex with deficit. Observe that the value of a max 
preflow equals the value of a max flow. We can limit the value of a max flow by some value d, 
by adding a new vertex s' and an arc (s',s) of capacity d and regarding s' as a source instead 
of s. 

An st-cut is a partition of the vertex set V into two sets S and V \ S, such that s S S and 
t G V \ S. The set of arcs going from S to V \ S is called the cut-set. The wa/ue of an st-cut is 
the sum of capacities of all arcs in the cut-set. The max-flow min-cut theorem |13] states that 
the value of the max st-flow is equal to the value of the min st-cut. 

Let M be some fixed value greater than the sum of all capacities in the graph, that is 
M > YleeE c ( e )- When we say that we assign infinite capacity to an arc, for example for vertex 
splitting, we in fact assign finite capacity of value M. This way we keep the flow assignment 
finite. 

2.2 Pieces and Recursive Decomposition 

Let G be a plane graph. We assume, for simplicity, that G is triangulated and has bounded 
degree. Both can be satisfied simultaneously using a standard vertex-splitting argument to 
get degree three and then triangulating each face such that each vertex degree in that face is 
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increased by at most 2. In Section [SJ where we report the cut-sets in time proportional to their 
sizes, we report the cut-set from the original graph without these simplification assumptions. 

A piece is a subgraph of G with the same embedding as G. A vertex of P is called a boundary 
vertex of P if it is incident in G to a vertex not belonging to P. All other vertices of P are 
called interior vertices of P. We let dP denote the set of boundary vertices of P. A hole of P 
is a face of P which is not a face of G. In certain places, we shall regard a hole as the subgraph 
of G embedded inside of it. It should be clear from context what is meant. 

A decomposition of a piece P is a set of sub-pieces Pi , . . . , Pj, such that the union of the 
vertex sets of these sub-pieces is the vertex set of P and such that every edge of P is contained 
in a unique sub-piece. We define every boundary vertex of P to be a boundary vertex of every 
sub- piece Pi that contains it. We change the standard definition of a decomposition a little and 
allow two sub-pieces to share edges. Edges that connect boundary vertices of a sub-piece are 
always included into this sub-piece, even if these edges belong also to another sub-piece. This 
does not increase the size complexity of the decomposition, and will simplify the discussion of 
the dual graph of the sub-piece. 

A recursive decomposition of G is obtained by first identifying a decomposition of G and 
then recursing on each sub-piece. The recursion stops when pieces of constant size are obtained. 
The recursive decomposition of G is the collection of pieces constructed over all levels of the 
recursion. We are interested in a special type of recursive decomposition satisfying the following: 
a piece P with p vertices and b boundary vertices is divided into a constant number of connected 
sub-pieces each containing at most \p vertices, at most ^b boundary vertices inherited from P, 
and at most 0{y/p) additional boundary vertices. In addition, we require that each piece has 
only a constant number of holes. When we refer to a recursive decomposition of G in the 
following, we shall assume that it is of this special type. Parent /child and ancestor /descendant 
relationships between pieces correspond to their relationships in the decomposition tree. 

Obtaining a recursive decomposition that satisfies all of the above conditions is non-trivial. 
The assumptions that subpieces contain at most \p vertices, at most ^b boundary vertices 
inherited from P, and at most 0(y/p) additional boundary vertices can be satisfied by finding 
an r-division of P, for r = ^p; see |21) for details. A construction of connected pieces, as we 
require, is given in [5]. In order to ensure that pieces are connected, we allow 0(n) new edges to 
be added while maintaining planarity and the chosen embedding of G. The arcs corresponding 
to these edges have zero capacity and thus do not affect max flows or min cuts. 

The results in the remainder of this subsection are taken from [5]. 

Lemma 1. A recursive decomposition as described above can be computed in 0(n log n) time. 

In order to bound the running time of our algorithms we need the following lemma. 

Lemma 2. Let V be the set of pieces in a recursive decomposition of G. Then ^2p e -p\P\ = 
0(n log n) and J2peV \dP\ 2 = O(nlogn). 

The proof of Lemma [2] follows by combining the following two lemmas with the fact that the 
depth of the recursive decomposition is O(logn). 

Lemma 3. Let Vi be the set of pieces in level i of a recursive decomposition of G. Then 
Epevj p \=0(n). 

Proof. Let c m \ n be the constant such that pieces of size at most c m \ a are not decomposed further 
in a recursive decomposition. We may assume that no piece has size less than |c m i n . Let L(p) 
denote the total number of vertices (counted with multiplicity) in the leaf-pieces of the recursive 
decomposition of an p- vertex piece P. We will show that L(p) < c\p—C2^/p for suitable constants 
ci and C2- The lemma will follow since the number of vertices (counting multiplicity) in any 
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level of the recursive decomposition is dominated by the number of vertices in the leaves; this is 
bounded by L(n) which is 0(n). 

We prove that L(p) < C\p — c 2y Jp by induction. In the base cases, in which a piece of size 
p is a leaf of the recursive decomposition (so p £ [^c m i n , c m i n ]), L(p) = p. This is bounded by 
cip - c 2y /p so long as c min < ci(±c min ) - c 2y Jc min . Setting c min > {c 2 /(\ci - l)) 2 guarantees 
this. 

Now consider a piece P of size p > c m i n . Assume inductively that the claim holds for all 
values smaller than p. Let pi be the size of the i th child of P; P has N = O(l) children. By the 
inductive hypothesis we get: 

l (p) = L ( p <) - ^2 ( ClPi ~ C2 VS) = ci ^Pi - c 2 (1) 

i i i i 

We lower bound ^2% ^/P~i by observing that ^/p. can only be as small as allowed by p < ^2iPi, 
Pi £ [0, §]. The minimum value occurs when two of the p^s are equal to | and all others are 
zero, giving: 

^vp,.>2wf=y^ (2) 

i 

Also note that when P is decomposed, only the new boundary vertices are replicated among the 
child pieces. Let us denote by a the constant of the planar separator theorem, then 

Y^\Pi\<\p\ + <rNy/^\. (3) 

i 

Combining Equations ([3j) , dH) , and (|2]), we get that 

L{p) < cxijp + (tNt/p) - c 2 y / 2p = cip - c 2 yV2 - ^aNj sjp 

This completes the induction for c 2 sufficiently larger than c\<jN . □ 

Lemma 4. Let Vi be the set of pieces in level i of a recursive decomposition of G. Then 
E P evJdP\ 2 = 0(n). 

Proof. Let P be any piece or the whole graph. By adding dummy boundary vertices that never 
become boundary nodes in the children, we may assume that: 

\dP\ > c^J\P\ (4) 

where c is a constant that we will pick below. In the case when P is the whole graph G we 
write dG do refer to these dummy vertices. We will show that for any piece P with children 
Pl,...,P N : 

X>Pif < |dif. (5) 

3 

The lemma follows from this because, by summing over all pieces in a level we get, 



£ \dP? < £ |^P| 2 <---< £ \dP\ 2 = \dG\ 
Pev t PeVi-i Pev 



2 



Since G has only dummy boundary vertices, \dG\ = c(y |G|) , which is 0(n), as desired. 
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We now prove Equation In the next equation, the first and second inequalities follow 
from the definition of recursive decomposition and Equation (J4j) , respectively: 

\dP j \<~\dP\ + ay/\P\<(~ + ~}\dP\ (6) 



Similarly, as in previous lemma we observe that only the new boundary vertices are shared 
between the children, so: 

\ dp i\ ^ \ dP \ + <rNy/\P\. (7) 

i 

Combining it all together we get: 

\dPj\ 2 < ( ~ + - J \dP\ Y \ dp j\ h Y Equation @ 

3 3 

< [ I + - J \dP\ (\dP\ + aNy/\P\) W Equation © 



£ ( - + - J \dP\ (\dP\ + criV-^ ) by Equation © 



< \dP\ 2 for constant c sufficiently large. 
This completes the proof. □ 



2.3 Nesting of Outer Boundaries 

Let V be the set of pieces in the recursive decomposition. We may assume that the fixed source 
s is not a boundary vertex of any piece of V. We can always add a new vertex s' to G with 
an infinite capacity edge to s and regard s' instead of s. This guarantees that the following 
definitions are unambiguous. If a hole H of P contains s, we refer to H as the outside of P. 
The subgraph of G contained in the outside of P is denoted ext(P). The outer boundary of P is 
the set of boundary vertices of P contained in the outside of P. Note that the outer boundary 
of a piece is contained in a single cycle of that piece. 

For simplicity, we also assume that each sink belongs to a piece having an outside. This is 
true for every vertex that belongs to a piece that s does not belong to, that is for all vertices of 
the graph except for a constant number of vertices which belong to the same leaf piece of the 
recursive decomposition as s. For this set of sinks we can compute max si-flow for every sink in 
0{n log n) total time using the algorithm of Borradaile and Klein [3]. 

These assumptions guarantee that the outer boundaries of the recursive decomposition has 
a rather simple structure as illustrated in Figure |2^a) . The outer boundaries nest and form a 
laminar family, i.e., a forest. Let the nesting of outer boundaries be represented by a forest P. 
The outer boundary C is an ancestor in P of an outer boundary C if the s-side of C contains 
the s-side of C. Each outer boundary C is represented in the forest only once, even if it is shared 
by multiple pieces. For an example, see Figure |2{b). 

For any C G T which is not a root in P, we define the outer piece of C as follows: 

1. If C bounds a hole H of a piece P £ V, which is not the outside of P, then the deepest 
piece P in the recursive decomposition which satisfies this condition is the outer piece of 
C (see piece a in Figure [3J . 
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Figure 2: Figure (a) shows a recursive decomposition that satisfies out all assumptions. In 
Figure (b), the corresponding forest J- is depicted. 




Figure 3: A single step of building the recursive decomposition. The piece is divided by finding 
an r-division. The outer piece of the outer boundary of a is b. The outer piece of the outer 
boundary of c is the part of a outside of c, which contains d, e, /, g, and the outer boundary of 
a. 



2. Otherwise, C does not bound a hole of any piece, except the outside of the pieces of which 
it is the outer boundary. Let P be a piece such that C is the outer boundary of P. Let C 
be the parent of C in T and let P' be the deepest piece in the recursive decomposition such 
that C is the outer boundary of P' . We define the outer piece of C to be the intersection 
of P' with ext(P) (see piece c in Figure [3]). Note that in this case the outer piece of C is 
not a piece of the recursive decomposition. 

To simplify our description, we shall assume in the following that all cycles C G J- have outer 
pieces of the first type. In particular, we assume that all outer pieces belong to the recursive 
decomposition. Dealing with outer pieces of the second type does not increase our time bound. 
We simply can include these additional outer pieces into V . The following lemma implies that 
handling such pieces does not affect the asymptotic running time of the algorithm. 

Lemma 5. The number of trees in {J^l is O(l) and each tree in J- has constant degree. 

Proof. By assumption, s is not a boundary vertex of any piece. For any outer boundary C 
corresponding to a root in J 7 , there is no other outer boundary in T separating C and s. Hence, 
every such C must lie inside a hole in the same piece P. Since P has only 0(1) holes, the first 
part of the lemma follows. 

For a node of J 7 to be a child of another node in J 7 , the outer boundaries corresponding to 
these nodes must share vertices with holes of the same piece. The second part of the lemma 
follows, again since a piece only has 0(1) holes and each piece is divided into a constant number 
of connected subpieces. □ 

2.4 Dense Distance Dual Graphs 

For a piece P, Fakcharoenphol and Rao |12] define the dense distance graph of P as the com- 
plete directed graph on the set of boundary vertices of P, representing shortest path distances 
among them in P. We apply a similar construction to the dual graph, guided by the recursive 
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decomposition. We define dense distance dual graphs for the pieces of the decomposition and 
for the holes of these pieces. We define a boundary face (of a piece or of a hole) to be a face 
incident to a boundary vertex (of the same piece or of the same hole). Our definition of dense 
distance dual graph is over the set of dual vertices corresponding to boundary faces. 

Let P be a piece in the recursive decomposition, denote by G*[P*] the subgraph of G* 
induced by the dual vertices corresponding to the faces of P. The internal dense distance dual 
graph IDDG*(P) of P is the complete directed graph over the set of dual vertices corresponding 
to boundary faces of P. An arc (/i, f2) in IDDG*(P) has weight equal to the (possibly infinite) 
shortest path distance from f\ to /2 in G*[P*]. Observe that due to our assumption that vertices 
have constant degree, the number of vertices in IDDG*(P) is 0(\dP\). The number of edges in 
IDDG*(P) is 0(\dP\ 2 ), and thus by Lemma [2j the total size of the internal dense distance dual 
graphs of all pieces of the recursive decomposition of G is 0(n log n). 

Let H be a hole of P. We define the dense distance dual graph of H, denoted DDG*(P), 
in a similar way. It is a complete directed graph on the set of dual vertices corresponding to 
boundary faces of H. Each arc (/i,/2) of DDG*(H) has weight equal to the shortest path 
distance in from dual vertex f\ to the dual vertex /2 in the subgraph of G* induced by the faces 
of H. Lemma [2] and our assumption that vertices have constant degree imply that the total 
number of edges of DDG* (H) over all holes of a single piece P is 0(\dP\ 2 ) and that the total 
number of edges over all dense distance dual graphs of all holes in the recursive decomposition 
of G is O(nlogn). We do not compute DDG*(P) for the hole H of P such that H = ext(P). 

The advantage of Fakcharoenphol and Rao's dense distance graphs is a faster implementation 
of Dijkstra's algorithm which we refer to as FR-Dijkstra. This algorithm can find a shortest path 
tree of a given vertex in a graph composed of dense distance graphs of pieces of a planar graph 
and additional arcs from the (original) planar graph in <3(61og 2 6) time, where b is the sum of 
the total number of vertices in the dense distance graphs, counted with multiplicity, and the 
total number of additional original arcs. We can use the implementation of FR-Dijkstra without 
any change also on dense distance dual graphs together with dual arcs from G*. 

Following Fakcharoenphol and Rao, we use FR-Dijkstra also for the construction of the 
dense distance dual graphs. Klein |27] gave a faster construction of dense distance graphs that 
construct dense distance graph in 0(|P| log |P|) time for a piece P with dP = 0(\/P), however 
we cannot use it for computing dense distance dual graphs of holes, since they need not obey 
this connection between the number of boundary faces and the size of the hole. We give more 
details about the algorithm of [27j below. 

Lemma 6. We can compute for all pieces P £ V the graph IDDG*(P) and the graphs DDG* (H) 
for every hole H of P other than ext(P), in total of 0(n log 3 n) time. 

Proof. First we compute IDDG*(P) for every PdP using Klein's construction [27| in 0(n log 2 n) 
total time. We proceed bottom-up on each tree T € T . Note that all the holes that we are 
interested in are bounded by some outer boundary in J 7 . For each C £ T, we consider the 
hole H inside C and compute DDG* (if) as follows. If C is a leaf of T then the construction is 
trivial, since the size of H is bounded by a constant. Otherwise, let Pc be a piece whose outer 
boundary is C. We build a graph D* composed of IDDG* (Pc) and dense dual distance graphs 
DDG*(H') for all holes H' ^ ext(Pc) of Pc- Note that for every such H', the outer bound- 
ary that bounds H' is inside the outer boundary C which bounds H. Thus we have already 
computed DDG*(H'). We compute DDG* (H) by running FR-Dijkstra in D* from each of the 
0(|C|) = 0(\dPc\) boundary faces of H. Each run takes 0(\dPc\ log 2 |9Pc|) time, so in total 
we require 0(\dPc\ 2 log 2 |<9Pc|) time for each outer boundary C. Applying Lemma [2] for the 
whole recursion we get 0(n log 3 n) total time. □ 

We will use FR-Dijkstra also for the rest of our algorithm. Using FR-Dijkstra with reduced 
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lengths defined by a potential function (see Section 14, ip will turn out to be very useful for our 
purposes. The original implementation of FR-Dijkstra by Fakcharoenphol and Rao |12| does not 
support reduced lengths. Therefore, we use an extension of FR-Dijkstra by Kaplan et al. |25j 
which allows using them. 

2.5 Multiple Source Shortest Path 

Another useful tool which we use in our paper is the multiple source shortest path (MSSP) data 
structure of Klein [27], which we already mentioned above. The MSSP data structure computes 
efficiently a set of distances such that all sources are different vertices along a single face of a 
planar graph. The data structure computes k distances in 0((n + k) logn) time. In addition, it 
can also produce an 0(n) space data structure, which reports the shortest paths themselves in 
0(1 log logn) time per path of length I. 

We use the MSSP data structure to compute dense distance dual graphs and for reporting, 
in Section [5j the dual paths that constitute parts of the cut-sets which we report. In a typical 
use of the MSSP data structure by our algorithm, we compute distances among boundary faces 
adjacent to a boundary cycle C of a piece P that bounds a hole He, inside the hole (we also use 
the MSSP for distances inside P; this case is similar). The internal part of the piece P is not 
part of Hq, and there is a single face in the embedding of He that corresponds to P. In the dual 
graph of He, all the dual vertices that we are interested in, which correspond to the boundary 
faces, are adjacent to the single vertex that represents the face of He that corresponds to P. 
Once we remove this vertex, together with the edges dual to C, which are all the edges adjacent 
to this vertex, all the vertices that are dual to the boundary faces are on a single dual face and 
we can use the MSSP algorithm on them. 

3 Single Source - All Sinks Max Flow 

For a fixed source s, we wish to find the value of a max si-flow in G for every sink t £ V \ {s}. 
In this section we describe an algorithm which is based on computing the value of a max flow 
from s to every outer boundary in J 7 in a recursive fashion. From this we can easily find the 
values we are looking for. In Section |4l we describe an efficient implementation of the algorithm. 

3.1 Max Preflow in a Separated Graph 

First, we consider a more specific problem, which enables us to implement the recursive step 
of our algorithm. Let G\ = (Vi,Ei) and G 2 = (V 2 ,E 2 ) be non-empty subgraphs of G, where 
Vx U V 2 = V, E 1 U E 2 = E and E\ n E 2 is a simple cycle C with vertex set V\ D V 2 . We shall 
identify C with V\ fl V 2 . Assume that s € V\\C and t G V 2 \ C. This assumption can be 
made without loss of generality. Algorithm MaxPreflowInSeparatedGraph finds a max 
st- preflow in G, given a max sC-preflow f s c in G\\ 

Algorithm MaxPreflowInSeparatedGraph(Gi , G 2 , f s c)'- 

1. Find a max Ct-flow fet in G 2 . Denote the value of this flow by det- 

2. Add the two flow assignments f s e and fet to form a pseudoflow in G. 

3. Send flow among vertices of C until there is no residual path from a vertex u £ C with 
excess to a vertex v G C with deficit; update the residual network accordingly. 

4. If there are vertices in C with remaining deficit, let —I be the sum of deficits over all 
vertices of C with a deficit; reset the flow on G 2 , and rerun the algorithm with a limit of 
det — ^ on the Ci-flow pushed in step [1] 
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The above algorithm is based on the flow partition scheme of Johnson and Venkatesan |24] , 
which was also used in Our implementation is a little different from the original scheme in 
its last step. First, we do not return excess flow from C to s; we are not required to do so, since 
we are looking for a preflow. Second, we do not return deficit flow from C to t. Instead we 
recompute a Ci-flow with value limited to dct — £■ This computation is equivalent, but as we 
show below, we can implement it faster. 

Notice that since edges of E\ n Pi are incident only to vertices of C, we may assume that 
both f s c and fct do not assign flow to any of these edges, since the vertices of C are all sinks 
of the first and sources of the second. Therefore at step [2] we indeed obtain a pseudoflow that 
does not violate the arc capacities. 

We shall refer to an algorithm implementing step[3l that is sending as much flow as possible 
from excess vertices of C to deficit vertices of C, as a flow fixing procedure along C in the residual 
network. 

3.2 Max Preflow in Pieces 

The main algorithm of this section, which we call Algorithm MAXPREFLOWToBoUNDARIES, 
uses Algorithm MaxPreflowInSeparatedGraph in recursion guided by the nesting of outer 
boundaries. The problem that Algorithm MAXPREFLOWToBoUNDARIES solves is to identify, 
for each piece P, a max preflow from s to the outer boundary of P. This max preflow is entirely 
outside of P. 

Instead of considering every piece, it is enough to consider only the unique outer boundaries 
in T, since if two pieces share the same outer boundary, we do not need to compute a max 
preflow from s to this outer boundary twice. Our algorithm traverses each tree T £ T from root 
to leaves. Let C be the outer boundary corresponding to the root node of T. The algorithm 
gets the max preflow f s Q from s to C outside of C as a parameter. This max preflow can be 
computed by embedding a super-sink t' inside C, connecting every vertex of C to t! by infinite 
capacity arcs, applying a max si'-flow algorithm for planar graphs, and removing t' . 

Algorithm MaxPreflowToBoundaries(T,C,/ s c): 

1. For every child C of C in T, perform the following steps: 

2. Add a super-sink t' inside C' connecting every vertex of C to t' using infinite capacity 
arcs. 

3. Apply Algorithm MaxPreflowInSeparatedGraph(Gi , G 2 , f s c) to find a max preflow 
f s f from s to t', where G\ is the part of the graph outside of C and G2 is the part inside 
C and outside C. 

4. Remove t' and the arcs incident to it, from f st i to obtain f s c', a max sC'-preflow. 

5. Apply MaxPreflowToBoundaries(T,C",/ sC -')- 

Algorithm MAXPREFLOWToBoUNDARIES considers every outer boundary C which is a 
child of C in T. The max preflow from s to C outside of C is given as a parameter of the 
procedure. The algorithm finds a max preflow from s to C", using a super-sink t' inside C 
and infinite capacity arcs as before. This way the problem becomes finding a max preflow from 
s to t'. (Note that C and C might share some vertices, in this case we do not connect the 
shared vertices to t', to avoid flow of infinite value.) We find the max preflow from s to t' using 
Algorithm MaxPreflowInSeparatedGraph. We can regard C as a separator in the graph, 
where G\ is the s-side (outside) of C and G<i is the i-side (where C' is). In the rest of this section 
and in the following section we focus on the fast implementation of the above steps. 
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Algorithm MAXPREFLOWToBoUNDARIES computes a max preflow from s to the outer 
boundary of P, in ext(P), for every piece P in the recursive decomposition. Every sink t is 
contained in some piece of constant size Pf. We find a max preflow from s to each sink t, using 
algorithm MaxPreflowInSeparatedGraph with the max flow f s c as a parameter, where C 
is the outer boundary of Pf. From this, we can obtain the value of a max si-flow for each sink 
t ^ s in G, as required, since this value is the same as the value of the max si-preflow that we 
found. 

3.3 Pushing Flow Through Holes 

Let C, C, and G' = G\ + G 2 be as defined in the description of algorithm MaxPreflowTo- 
BOUNDARIES. Let P be the outer piece of C . Note that C is the outer boundary of P. In 
step CD of Algorithm MaxPreflowInSeparatedGraph we push flow from C to C in G 2 . A 
simple way to implement this step is using the max flow algorithm of Borradaile and Klein [3] in 
Gi with a super-source connected to C and a super-sink connected to C'\ C, in 0(| G2 1 log |G 2 |) 
time. If P has no holes except of the holes bounded by C and by C', then G2 is P, and we 
can use the simple solution. However, if there are more holes in P, other than these two, then 
G2 might be much bigger than P, since it should contain also the parts of the graph inside the 
additional holes. We cannot afford to spend that much time in this step. Rather, we would like 
to bound the amount of work here by the size of P, so that we can use Lemma [2] to bound the 
total time of our algorithm. 

Let Hi, . . . , Hk denote the holes of P bounded by neither C nor C and let Ci, . . . , be 
the boundaries of these holes. We observe that G 2 = P + H± + • • • + H^. In this subsection, we 
show how to decompose the computation of the max flow found in step [1] of Algorithm MAX- 
PreflowInSeparatedGraph into smaller steps according to this partition of G2, such that 
these steps consist of max flow computations inside P and of calls to the flow fixing procedure. 
In Section U we will give an efficient implementation for these steps. 

The single additional hole case We first consider the case where k = 1. Let H = H\ 

and let Ch = Ci be the boundary of H. We give here an algorithm PushFlowThrough- 
Holes(P, C, C , H) which finds a max flow inside P + H from C to C'. The algorithm is as 
follows (each step uses the residual network for the flow found in earlier steps): 
Algorithm PushFlowThroughHoles(P, C, C ,H): 

1. Compute a max flow in P from C to C. 

2. Compute a max flow in P from C to Ch- Denote the value of this flow by d 2 . 

3. Compute a max flow in P from Ch to C . Denote the value of this flow by d^. 

4. Run a flow fixing procedure along Ch in P + H. 

5. Let £2 be the sum of excesses over vertices of Ch with an excess, and let —£3 be the sum 
of deficits over all vertices of Ch with a deficit; reset the flow to zero in P + H, and rerun 
the previous steps with a limits of d 2 — ^2 an d ^3 — £3 on the flow pushed in steps 2 and 
3, respectively. 

Lemma 7. Algorithm PuSHFLOWTHROUGHHOLES(P, C, C, H) finds a max flow from C to 
C inP + H. 

Proof. We will show that after (the first execution of) step 4, we have a pseudoflow with no 
residual paths in P + H from C to C", from C to deficit vertices on Ch, from excess vertices on 
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Ch to C, or from excess vertices on Ch to deficit vertices on Ch- This suffices to show that at 
termination we have a max flow from C to C in P + H, since step 5 in fact only returns excess 
to s and deficit to t, as in the flow partition scheme of Johnson and Venkatesan [24) . 

Let us identify a cut with the corresponding cut-set, that is the set of arcs crossing the cut. 
After step 1, there is a saturated cut K\ in P separating C from C. Thus, any residual paths 
from C to C' in P + H must cross H . After step 2, such residual paths cannot exist. Note that 
K\ stays saturated after step 2. We now have a new saturated cut K2 in P with C on one side 
and C and Ch on the other side. After step 3, K2 stays saturated (since C and Ch are on the 
same side of K<±) and we get another saturated cut K3 in P separating Ch and C. 

Consider an augmenting path implementation of the flow fixing procedure in step 4. Each 
augmenting path has both its endpoints on Ch so it crosses neither K2 nor K3 (as it would 
have to cross K2 or K3 in both directions). Hence, K2 and K3 stay saturated after step 4 so 
there is no residual path in P + H from C to C, from C to Ch, or from Ch to C". The flow 
fixing procedure ensures that there are no residual paths in P + H from excess vertices on Ch 
to deficit vertices on Ch- This shows the desired. □ 

Generalizing PushFlowThroughHoles to k additional holes Now let us generalize 
Algorithm PushFlowThroughHoles(P, C, C", H) to arbitrary k. Instead of a single hole 
H, it gets a set {Hi, . . . , H k } of holes as its fourth parameter. For k > 2, we can regard 
P + Hi + . . . + -ff/c-i as a piece P' with one hole H k in addition to the two bounded by C and 
C' . Hence the call PushFlowThroughHoles(P, C, C , {Hi, . . . , H k }) solves the problem. 
Plugging in these parameters in the procedure for the single hole case, we get the following 
recursive algorithm: 

Algorithm PushFlowThroughHoles(P, C, C',{H 1} H k })\ 

1. Execute PushFlowThroughHoles(P, C, C, {Hi, . . . , Hf.-i}) (finds max flow in P' 
from C to C). 

2. Execute PushFlowThroughHoles(P, C, C k , {Hi, . . . , H k _{\). Denote the value of this 
flow by d,2 (finds max flow in P' from C to C k ). 

3. Execute PushFlowThroughHoles^, C k , C, {Hi, . . . , iJ fe _ 1 }). Denote the value of 
this flow by (I3 (finds max flow in P' from C k to C). 

4. Run a flow fixing procedure along C k in P + Hi + . . . + H k (i.e., in P' + H k ). 

5. Let £2 be the sum of excesses over vertices of C k with an excess, and let —I3 be the sum of 
deficits over all vertices of C k with a deficit; reset the flow to zero in P' + H k , and rerun 
the previous steps with a limits of d>2 — I2 and (I3 — £3 on the flow pushed in steps 2 and 
3, respectively. 

Lemma 8. A call to PushFlowThroughHoles(P, C, C , {Hi, . . . ,H k }) computes a max 
flow from C to C in P + Hi + . . . + H k . The total number max flow computations in P between 
boundaries of holes and the total number of calls to a flow fixing procedure is 0(1) (for constant 
k). 

Proof. Correctness follows from Lemma [7] using induction on the recursive procedure. The 
number of max flow computations in P and the number of calls to the flow fixing procedure are 
both exponential in k. Since k = 0(1), the second part of the lemma follows. □ 
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4 An Efficient Implementation 



In this section, we give an efficient implementation of Algorithm MaxPreflowToBound- 
ARIES. Recall that this algorithm applies Algorithm MaxPreflowInSeparatedGraph to 
obtain a max preflow from s to an outer boundary C', outside of C' , given a preflow from s to 
C, outside of C, where C is the parent of C in T . We show that by maintaining flow assignments 
implicitly, such a call to Algorithm MaxPreflowInSeparatedGraph can be implemented 
to run in 0((\P\ + |<9P| 2 )log 2 \P\) time, where P is the outer piece of C . Combined with 
Lemma |21 this gives us the required 0(nlog 3 n) time bound for Algorithm MaxPreflowTo- 
BOUNDARIES. An important tool that we use is the fast implementation of the flow fixing 
procedure by Borradaile et al. |3] which we go through briefly in the following. 

4.1 The Flow Fixing Procedure of Borradaile et al. 

Recall that the flow fixing procedure gets as input a cycle separator C for a graph G = G\ + G2 
and a pseudoflow in G where all excess and deficit vertices (that are not at sources or at sinks) 
are on C. The procedure updates the flow network by sending as much flow as possible, among 
vertices of C, from excess vertices to deficit vertices. The efficient implementation of the flow 
fixing procedure from [3] takes advantage of the fact that the vertices of C are on a single simple 
cycle in the plane. It processes the vertices of C one by one in cyclic order. At each step, if 
the current vertex Vi has excess, as much flow as possible is sent from Vi to the yet unprocessed 
vertices. Similarly, if V{ has a deficit, as much flow as possible is sent from the unprocessed 
vertices to V{. By processing the vertices one by one, the procedure maintains the invariant that 
no flow can be sent from a processed vertex of C with excess to any other vertex of C, or to a 
processed vertex of C with deficit from any other vertex of C, without violating the capacity 
constraint of the network. 

An important tool for the fixing procedure is an algorithm of Hassin [19j for (single-source 
single-sink) max flow in a plane digraph where the source and sink are on the same face. This 
algorithm works by adding an infinite capacity arc eoo from the sink to the source and then 
computing shortest path distances in the dual graph from the face to the right of eoo. These 
distances define a potential function over the face which induces a max flow in the primal graph. 
Specifically, the flow /(e) on an arc e is defined by the difference between the potential of the 
face to the left of e and the face to the right of e. This flow assignment defines a circulation, 
and after removing 600 from this circulation we get a max st-flow. 

In the flow fixing procedure we connect the unprocessed vertices of C with edges of infinite 
capacity in both direction, thereby essentially identifying all these vertices with fj+i, the vertex 
that follows Vi on C. Then we can use Hassin's algorithm to compute a max flow from Vi to 
if Vi has an excess, or from to Vi if Vi has a deficit. This way using 0(|C|) calls to Hassin's 
algorithm we send flow among the vertices of C such that at termination there is no residual 
path from any excess vertex C to any deficit vertex of C. 

The fixing procedure uses infinite capacity arcs between consecutive vertices of C. There is an 
arc eoo between Vi and Uj+i as well as infinite capacity arcs between two consecutive unprocessed 
vertices of C. The vertices of C are connected by a cycle anyway, and we do not want to change 
the graph G during the fixing procedure. Therefore, instead we change the capacities of the 
relevant arcs of C to infinity when needed. Instead of removing arcs with infinity capacity, we 
change back the capacity of the arc to its previous value. Let e be an arc of C such that we 
changed the capacity of e to infinity and back to its previous value. The potential function that 
we found using Hassin's algorithm defines a flow of some value on e. If e is the arc between Vi 
and then we should remove the flow assigned to e in order to get the correct max UjVj + i-flow 
(or max Vi+iv j-flow). In addition, since we should send as much flow as possible from Vi to Vi+\ 
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if Vi has an excess, or in the opposite direction of v- l+ \ has a deficit, we assign to rev{e) the 
maximum amount of flow that its residual capacity permits. If e is one of the arcs that connect 
unprocessed vertices of C we should remove the flow assigned to e, so we will not violate its 
capacity. This way, the result of each iteration of the fixing procedure is a potential function for 
the faces of G and a flow assignment for the arcs of C which we should add or reduce from the 
circulation induced by the potential function. 

During the course of the fixing procedure, we maintain the current flow among the vertices 
of C by accumulating the flow assignment for arcs of C and a potential function (j) over the face 
of G, which defines a circulation on this graph. Every iteration uses the residual network for 
the flow found in all previous steps. 

As we described, Hassin's algorithm finds a max flow by computing a shortest path tree in 
the dual. To describe this part further, let us define X* to be the set of faces incident to vertices 
of C. This set contains the dual vertices from which we start the shortest path computations. 
Since the vertices of G have constant degree, we have that \X*\ = 0(\C\). Moreover, let G*_ c be 
the graph obtained from G* by removing the dual arcs of C and their reverses. Note that this 
splits G*_q into two disconnected parts, the inside and outside of C. Each iteration of the fixing 
procedure applies Hassin's algorithm, which compute shortest path distances from a of X* in 
G* with respect to a weight function on the arcs induced by the arc lengths l(u, v) of G* and the 
current potential function (p. More precisely, the reduced length 1^ with respect to <j) is used, where 
l^{u,v) = l(u,v) + 4>{u) — 4>(v) |23j . A crucial observation to make in order to get an efficient 
implementation is that the only reduced distances that we are interested in in G* are those 
starting and ending in vertices of X* and these can be obtained from the X*-to-X* distances in 
G*_ c and the restriction of (f) to X* , together with the lengths of the arcs dual to arcs of C. To 
see this, consider a path Q = V\Vi . . . v k in G*_ c with vi,Vk G X* . Then by a telescoping sums 

argument, the reduced length l<f,(Q) of Q is J2i=i K v i v i+i)-<ft{vi+i)+(f>{ v i) = KQ)+4>{vi)-4>{vk)-, 
showing the desired. 

In order to compute shortest path distances in G* , FR-Dijkstra is applied, where G*_q is 
represented by a matrix of X*-to-X* distances and the dual arcs of C are given explicitly. The 
algorithm only computes distances among vertices belonging to X* . As we saw above, these 
values suffice both to update the flow on C and to update the X*-to-X* reduced lengths in 
G*_ c . Hence, it suffices for the algorithm to maintain only the restriction of the face potential 
to vertices of X* . 

After processing all vertices of C, we get a restriction of the potential function (j) which 
defines the reduced lengths among vertices of X*. We want to extend this potential function 
to all faces of G, so that it will define a circulation in G. This circulation, together with the 
flow on arcs of C give us the desired flow among the vertices of C. Borradaile et al. |1] showed 
that using one call to FR-Dijkstra on X*, followed by one call to Dijkstra's algorithm over the 
entire graph G, we can get the desired potential function for all faces of G. Once we have this 
function, it is easy to compute the circulation it induces and add the flow on arcs of C to get 
the desired flow which balance the excesses and deficits at vertices of C as much as possible. 

4.2 Implementation of Algorithm MaxPreflowInSeparatedGraph 

We are now ready to describe our implementation of Algorithm MaxPreflowInSeparated- 
GRAPH. Let P, C, and C' be defined as in the beginning of the section and assume that a 
max preflow f s c from s to C has been computed. Denote by He and Hq> the outside of C 
and C', respectively. The algorithm gets the max sC-preflow outside of C, f s c, as an input 
and combines this preflow it with a max CC'- flow inside C and outside of C', to get a max 
sC'-preflow. We are not interested in the assignment of flow outside of C. In fact, we cannot 
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even represent this flow if we want to implement the algorithm withing the desired time bound. 
We only need a representation that allows us to run the efficient implementation of the fixing 
procedure described above. From the description of the fixing procedure above, we conclude 
that it is enough to have an implicit representation of the input preflow f s e which consists of 
two parts: 

1. The distances among the faces adjacent to C in ext(P), that is ~D~DG*(Hc), where the 
length of a dual edge is the residual capacity of the primal edge with respect to f s e- 

2. The excess at each vertex of C with respect to f s e- 

This does not cover the assignment of flow to arcs of C. However all vertices of C are sinks of 
f s C, therefore we may assume that the arcs of C do not carry any flow in f s e- This information 
is enough to run the efficient implementation of the fixing procedure described above. 

The output of Algorithm MaxPreflowInSeparatedGraph is a max sC'-preflow. If C is 
a single vertex (that is, a sink t), then we are only interested in value of the preflow. Otherwise, 
we use this output as an input for computing a max sC"-preflow, where C" is a child of C in P. 
Therefore, it is enough to use an implicit representation of this kind also for the max sC'-preflow 
that we find as an output. 

For an outer boundary C that is a root of a tree in P we need to produce an implicit 
representation of f s e from an explicit max sC-flow. We first find such a flow using the algorithm 
of Borradaile and Klein |3], and then we compute the dense distance dual graph DDG* (He) 
with respect to the residual capacity using the MSSP algorithm of Klein [27J. This is done in 
0(n log 2 n) time for all trees of P. 

Our goal is to implement Algorithm MaxPreflowInSeparatedGraph to run in 0((\P\ + 
|(9P| 2 )log 2 |P|) time. Let us first consider the simple case where the only holes of P are He 
and the hole bounded by C". This is the case where we do not require Algorithm PUSH- 
FlowThroughHoles from Section [3J3] In this case, Algorithm MaxPreflowInSeparat- 
EDGRAPH simply sends a max flow from C to C in P and then runs the flow fixing procedure 
on C outside C. 

The first part, sending a max flow fee from C to C in P, can be done in 0(|P| log \P\) time 
by a single call to the algorithm of Borradaile and Klein j3] in P (after adding a super-source and 
a super-sink, again we do not connect the shared vertices of C and C, if any, to the super-sink). 
From the max CC'-flow we computed in P we obtain IDDG*(P), the internal dense distance 
dual graph with distances representing the residual capacity. Using Klein's MSSP algorithm |27) 
it takes 0((\P\ + \dP\ 2 ) log |P|) time to construct IDDG*(P). 

We now have dense distance dual graph on DDG* (He) that represents f s c, and the dense 
distance graph IDDG*(P) that represents fed together with the amount of excess (from f s e) 
and deficit (from fee) a t each vertex of C. These dense distance dual graphs, together with 
the dual arcs of C, allow us to apply the flow fixing procedure of Borradaile et al. [1] to C in 
the part the graph outside C . However, in the last step of the flow fixing procedure, where the 
procedure extends the potential function defined over the face adjacent to C to all faces outside 
C, using a call to FR-Dijkstra and a call to Dijkstra's algorithm, we limit this extension only 
to P. On one hand, we cannot afford extending the flow to the entire subgraph outside C since 
it might be much bigger than P, on the other hand we do require to extend the flow to P, since 
we are interested in the flow that enters C which is not represented explicitly in the current 
representation. So instead of using Dijkstra's algorithm to extend the flow to the entire graph 
as is [3], we use another call to FR-Dijkstra, where the arcs dual to P and the arcs dual to C 
are represented explicitly, but He is represented using DDG* (He)- This allows us to run the 
flow fixing procedure and to extend the flow to P in 0((|P| + |<9P| 2 ) log 2 |P|) time. 
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If there are no deficit vertices on C, we have the desired s-to-C" preflow, represented explicitly 
in P and implicitly outside C. Otherwise, we rerun the above steps with a limit on the amount 
of flow pushed in P as described in step g] of Algorithm MaxPreflowInSeparatedGraph. 
In fact, in this case we can skip the computation of the explicit representation of the flow in P 
following the first call to the flow fixing procedure. 

Finally, we obtain an implicit representation of the combined flow f s c> using the MSSP 
algorithm. The total time for all these steps is 0((|P| + |<9-P| 2 ) log 2 \P\), as required. 

We assumed above that P has only two holes, Hq and the hole bounded by C. We now show 
how to handle the general case with k additional holes Hi, . . . , H^ bounded by cycles C±, . . . , Cfe, 
respectively. We execute step Q] of Algorithm MaxPreflowInSeparatedGraph using an 
efficient implementation of Algorithm PushFlowThroughHoles(P, C, C , {Hi, ... , H k }). It 
follows from the description of this algorithm and from Lemma [SJ that we should support a 
sequence of 0(1) operations of the following two types: 

1. Find a max flow in P from a cycle d to a cycle Cj (the single additional hole case of 
Algorithm PushFlowThroughHoles), 

2. Run the flow fixing procedure along Ci in P + Hi + . . . + Hi, for some 1 < i < k. 

During the course of Algorithm PUSHFLOWTHROUGHHOLES, we represent the current pseud- 
oflow in a way similar to the one above. We represent the pseudoflow explicitly in P, including 
the arcs of Ci, . . . , C}-, and implicitly for the arcs inside the holes Hi, . . . , H^, using a dense dis- 
tance dual graph for these holes with respect to the current residual capacity and the balance of 
flow at vertices of Ci, . . . , C k - Before the initial call to Algorithm PUSHFLOWTHROUGHHOLES, 
the flow in Hi, . . . , H^ is all zero, so we can represent the initial residual network in each Hi 
implicitly using DDG* (Hi) with the distances defined by the original capacity function. By 
Lemma [6j all of these dense distance dual graphs can be precomputed within the desired overall 
time bound of our algorithm. 

The first kind of operation, sending max flow from one cycle to another inside P can be 
implement by simply using the algorithm of Borradaile and Klein [3] with a super-source and a 
super sink in 0(|-P| log |P|) time. 

The second type of operation, running the fixing along d in P + Hi + . . . + iTj is done in 
a way similar to before, but this time we should also take into account the part of the graph in 
the other holes Hi, . . . , -ffj—i. Again, to apply the fixing procedure to Ci we require the dual 
arcs of Ci together with DDG* (Hi) and IDDG*(P). This time, in addition we also need the 
dual arcs of Cj for every 1 < j < k and also DDG* (Hi) for every 1 < j < i — 1. All of these 
dense distance dual graphs, as well as the length of the dual arcs of Cj's, are with respect to 
the current residual capacity. When we apply FR-Dijkstra during the fixing procedure, we do 
it over a graph that is defined by the union of all of these dense distance dual graphs and dual 
arcs of cycle separators. The total number of vertices remains 0(\dP\), so the running time of 
the flow fixing procedure remains 0((\P\ + |aP| 2 )log 2 \P\) as before. 

After each time we apply the fixing procedure, we update the residual capacity of the arcs 
dual to each cycle separator Cj and keep the restriction of the potential function for each of 
the dense distance dual graph DDG(Hj), so that it will represent correctly the current residual 
capacity. 

When we are done, we get the explicit flow of P from the dense distance dual graphs and 
the flow assignment on the arcs of the cycles, and from them we get the implicit representation 
of the max Ci-flow in the same way as we described before, without computing the explicit flow 
for arcs that are not in P. 

We have shown that the two types of operations required for finding the max Ct-flow take 
0((|P| + |<9-P| 2 ) log 2 \P\) time. By Lemma [81 there are only 0(1) of them during a call to 
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Algorithm PUSHFLOWTHROUGHHOLES, and hence step [T] of Algorithm MaxPreflowInSep- 
ARATEDGRAPH runs in 0((|P| + |<9P| 2 ) log 2 |P|) time. The fixing procedures at step[3]and step 
[4] are implemented in a similar way, within the same time bound. 

We can now conclude this section with the main result of the paper. 

Theorem 1. Given a planar n-vertex digraph G = (V,E) with a fixed source s in G, the max 
st-flow values from s to every sink t G V \ {s} can be computed in a total of 0(n log 3 n) time. 

5 Reporting Min Cut Sets 

The same ideas we used in Section [3] and Section [4] can also lead to an efficient algorithm for 
reporting minimum si-cut sets for a fixed s and a given t. For simplicity we shall refer to cut-sets 
simply as cuts. The algorithm is based on a similar preprocessing scheme as the algorithm from 
sections and [|] We use the duality between cut-sets in a planar graph and cycles in the dual 
planar graph. Let K be an si-cut set, that is a set of arcs whose removal separates s from t. 
Then, the dual arcs of the arcs of K form a cycle which separates s from t and goes clockwise 
around t. In particular, if / is a max st-flow, then a cycle M which separates s from t, and 
goes clockwise around t, whose length is equal to the value of / (which we denote by |/|) is a 
dual of a min-cut. In this case, we call M a shortest st-separating cycle. If we consider the dual 
of the residual graph (the lengths are given by residual capacities with respect to /), then the 
length of M is 0, since it corresponds to a saturated si-cut. The following observation is the 
key for this step. Fix a source s, a cycle C, and two faces x, y adjacent to C. Then, for every 
shortest st-separating cycle M (for any sink t), which contains a subpath Q from x to y that is 
embedded outside C, the path Q is one of two possible paths: either the shortest x to y path 
that goes clockwise around s, or the shortest x to y path that goes counterclockwise around s, 
regardless of the choice of t. Our algorithm, after some preprocessing, is capable of reporting 
the arcs of the shortest st separating cycle in time that is roughly proportional to its length. 

We begin by giving a simple characterization of si-cuts. For this characterization we assume 
that for every sink t, the value of the min si-cut is positive. It is easy to identify sinks that satisfy 
this condition, since these are the vertices reachable from s after removing all zero-capacity arcs. 

Lemma 9. Let f be an st-flow in G. Let M be a cycle of length in G*j, the dual of the residual 
network for f . Then, if there is any flow crossing M , or in other words if there is an arc e of 
G such that the dual arc of e is in M and /(e) ^ 0, then M is a shortest st-separating cycle. 
On the other hand, a shortest st-separating cycle M has length in G*j and there is some flow 
crossing M. Furthermore, for every arc e of G whose dual in M, /(e) > 0. 

Proof. Let M be a cycle of length in G*j, M separates the vertex set of G into two sets - S 
and T. Let S and T be such that every e = (u, v) whose dual is in M satisfies u G S and v G T. 
Since the total length of M is zero, we get that the residual capacity of every arc e whose dual 
is in M is zero. Therefore, for every such e, we get /(e) > 0. If there is flow crossing M, then 
the net amount of flow going from S to T is positive. Hence, it has to hold that s G S and 
i G T, which means that M corresponds to an si-cut. The residual capacity of M is zero, so it 
is a shortest si-separating cycle. 

The other direction follows from the discussion above, since a shortest si-separating cycle 
corresponds to a saturated si-cut. □ 

We use an algorithm similar to the one from the previous sections. However, this time finding 
a maximum preflow is not enough, since a max st-preflow might have a saturated cut whose 
value is greater than the value of the max si-flow. Therefore, this time we implicitly compute a 
max si-flow. 
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In contrary to the previous sections, where we guided our flow partition using a recursive 
decomposition, this time we use an r-division of G as defined in |15] , It is a decomposition of G 
into O(-) pieces of size 0(r). Each piece has 0{y/r) boundary vertices and a constant number 
of holes. An r-division can be computed in 0{n log r + {n/^/r) log re) time [21] . 

For every vertex t in the same piece as s, we compute a max si-flow and find a min st-cut 
using the algorithm of Borradaile and Klein [3], for a total of 0(rn log re) time. We store all 
these cuts explicitly using 0(rn) space. It is left to consider only the case where s and t are in 
different pieces. 

First, for every piece P in the r-division with outer boundary C, we compute a max sC-flow 
using the algorithm of Borradaile and Klein |3] (using a super-sink). We will use these max 
flows for our flow partition algorithm. The total time for these computations over all pieces is 
0(n 2 logn/r), and we produce implicit representations of these flows in total of 0(re 2 logre/r) 
time using the MSSP algorithm of Klein [27J. We also compute the dense distance dual graph 
of all holes in the r-division. Each of them is computed with the MSSP algorithm in 0{n\ogn) 
time, which results in 0(n 2 log n/r) time for all holes. 

Then, we apply Algorithm IMPLICITMAXFLOW which is described next for every sink t that 
is not in the same piece s. We denote by P the piece containing t, and let C be the outer 
boundary of P. This algorithm finds a max flow f s t from s to t, the flow is represented explicitly 
in P and implicitly in the rest of the graph. 

Algorithm iMPLldTMAxFLOW(i): 

1. Let fi be the max sC-flow strictly outside of C computed in advance (represented implic- 
itly). 

2. Compute max Ct-flow fi inside of C. The flow fi is represented explicitly in P, if there 
are additional holes inside C, then fi is computed using Algorithm PushFlowThrough- 
HOLES and fi is represented only implicitly in these holes. 

3. Combine f\ and fi, and run the flow fixing procedure along C. 

4. Compute the value of the max si-flow, denoted by d as follows. Let —I be the sum of 
deficits over all vertices of C with deficit. Then d := I/2I — L 

5. Restore the flow f\ from step 1. 

6. Push |/i I — d units of flow from vertices of C with excess back to s in the residual network of 
fi (strictly outside of C), thus obtaining an sC flow f[ of value d (represented implicitly). 

7. Compute a max Ct-flow inside C with a limit of d units on its value. 

8. Combine the flows f[ and f 2 and run the flow fixing procedure along C into a single flow 
fst, represented explicitly in P and implicitly in the rest of the graph. 

For every t, step 2 takes 0(r log 2 re) time using the max flow algorithm of Borradaile and 
Klein and Algorithm PushFlowThroughHoles. Step 7 is similar. In step 3 and step 8 we 
run the flow fixing procedure, which requires 0(r log 2 n) as well. It remains only to bound the 
time for step 6. 

In step 6, we return excess flow from C to s, which is a standard step in the flow partition 
scheme that we use. Johnson and Venkatesan |24j implemented the flow return step explicitly, 
we cannot do this since it will take too much time. Moreover, in our case the flow f\ is repre- 
sented implicitly. The solution we used in Algorithm MaxPreflowInSeparatedGraph was 
to recompute the flow, in our case /1, however we cannot afford spending the time required by 
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this recomputation here. Note that we compute f\ only once for all vertices of P. In contrast, 
in step 7 we send back flow from C to t as in Algorithm MaxPreflowInSeparatedGraph, 
we can do it in this step since recomputing a flow from t to C is possible within our time 
bound. Next we describe how to push back flow from C to s in the implicit representation of 
fi, obtaining an implicit representation of /{. 

It makes the implementation of this step easier if we assume that every face that is adjacent 
to a vertex of C is also adjacent to an edge of C. This assumption can be satisfied by vertex- 
splitting along all boundaries of the r-division, as we described in Section \2. 21 to obtain a graph 
of constant degree. After the vertex-splitting, the graph G is no longer triangulated, and the 
size of the boundaries might increase. However, the blow-up in the size of the boundaries is by 
a constant factor, since the degree of each vertex was constant also before the splitting. For a 
face x adjacent to C we denote by e(x) the edge of C that is adjacent to x. 

We divide the — d units of flow, which we want to push back from C, among the vertices 
of C in an arbitrary way, such that we will not return from any vertex amount of flow greater 
than the excess it has in f\. 

We show how the implicit representation of the flow outside C changes, when we push flow 
from vertices of C back to s, in the residual network of f\. Since the flow f± goes from s 
to C, the flow is pushed back along paths that do not contain any vertex of C, except their 
starting vertices. Recall that the implicit representation contains all pair distances between faces 
adjacent to vertices of C in the dual of the residual graph for f\. We first show the change in 
this implicit representation after pushing back k units of flow from some vertex v of C to s. As 
we show below, we get the same implicit representation regardless of the exact path R v along 
which the flow is pushed back. 

As in Section |4j we denote by Hq the hole of P outside C and DDG* (Hq) is the dense dis- 
tance graph representing the residual capacity in this hole. Consider an arc (x, y) in DDG*(-ffcO 
that connects two faces adjacent to vertices of C. It corresponds to a dual path Q in the dual 
of ext(-P), embedded strictly outside C. Assume that we push k units of flow along a path R v 
from a vertex v £ C to s. 

We say that a directed dual path M in the dual graph crosses a directed path N in the 
primal graph if M contains a dual of an arc of N or its reverse. In the former case, we say that 
M crosses N from the right side, whereas in the latter we say that M crosses N from the left 
side. If Q crosses R v exactly once, then after pushing k units of flow along R v - 

1. if Q crosses R v from the right side, that is it contains a dual of an arc of R v , then the 
length of Q decreases by k, 

2. if Q crosses R v from the left side, that is it contains a reverse of a dual of an arc of R v , 
then the length of Q increases by k. 

This can be further generalized. Consider the order of the crossings between Q and R v as we 
encounter them when we traverse R v . Consecutive crossings alternately increase and decrease 
the length of Q by k. If the total number of crossing is odd, then the length of Q changes by 
+k or —k, depending on the direction of the first crossing. On the other hand, if the number of 
crossings is even, then the length of Q remains unchanged (see Figure |4]). 

Observe that Q splits ext(P) into two regions E\ and E%. Let E\ be the region containing 
s. Formally speaking, Q does not start at C, but at two faces that are adjacent to vertices of 
C. We divide the graph outside of C with a curve that is composed of the embedding of Q, as 
well as two curves that connects x and y to e(x) and e(y) respectively. In particular, the cycle 
C itself is divided into two parts C\ and C2, where Cj belongs to E^. 

In the beginning of our algorithm, we computed DDG* (He) that represents the flow f\ 
outside C. We now extend this computation to make it possible to update the length of Q. Fix 
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(a) Q crosses R v an even 
number of times. 



(b) Q crosses R v an odd 
number of times. 



Figure 4: Illustration of lemmas [TT1 and [T2l The cycle C is marked in bold, Q is solid and R v is 
marked with a dashed line. If k units of flow are pushed along R v , then every crossing between 
Q and R v contributes a change of +k or —k, alternately, in the length of Q. 

an arbitrary vertex w on C an let R' be an arbitrary the path connecting w to s outside of C. 
Instead of computing DDG* (He), we compute for each pair of faces x, y that are adjacent to 
vertices of C and lie outside C two shortest paths connecting them. One of them is the shortest 
dual path that crosses R' an even number of times, while the other is that shortest dual path 
that crosses R' an odd number of times. It is clear that the shortest path from x to y, whose 
length defines the length of (x, y) in DDG* (He), is always one of these two dual paths. We show 
that the lengths of such x-to-y dual paths can be efficiently updated after some flow is pushed 
via paths from C to s. 

Lemma 10. Let Q be a dual path from x to y and let R v be a path connecting v S C and s, 
both strictly outside C . 

If v is in the same region as w, the vertex which defines R' , then the parity of the number 
of crossings of Q with R' is the same as the parity of the number of crossings of Q with R v . 
Otherwise, if v and w are in different regions, then these parities are different. 

Proof. We claim that the parity of the number of crosses of Q with R v depends solely on whether 
v is in C\ or Ci. Indeed, in every time R v crosses Q, it switches from E\ to E2 or vice versa. 
From this, we get the desired. □ 

From this lemma we infer that for every v E C\ the number of crossing between the path R v 
and the dual path Q is even, and for v G C2 the number of crossing is odd. 

Lemma 11. Let Q be a dual path from x to y strictly outside C . Assume that we push k units 
of flow from v E C to s along a path R v . Assume also that Q crosses R v an even number of 
times (see Figure \^(a)ty . Then, pushing the flow along R v does not affect the length of Q. 

Proof. Every time R v crosses Q, it switches from E\ to Ei or vice versa. Hence, Q alternates 
between dual arcs of R v and reverse of dual arcs of R v . A crossing in one direction decreases 
the length of Q by k while a crossing in the other direction increases the length of Q by k. The 
number of crossing of each type is the same, so the length of Q remains intact. □ 

The following lemma considers the clockwise order of x, v and y around C. The faces x and 
y do not belong to C, but their order around this cycle, relative to v, is well-defined by e(x) and 

e(y)- 



20 



Lemma 12. Consider the scenario from the previous lemma, but assume that Q crosses R v an 
odd number of times (see Figure ^b)^ . If x, v and y are in clockwise order around C then after 
pushing the flow along R v , the length of Q increases by k. Otherwise, if the order of x, v and y 
is counterclockwise, the length of Q decreases by k. 

Proof. We consider the case, when x, v and y lie in clockwise order around C. The other case 
is symmetric. 

In this case the crossing of Q and R v which is closest to v in R v is from the left side of R v . 
The following crossings alternate between crossings from the right side and crossing from the 
left side. This implies that Q crosses R v from the left side t + 1 times and from the right side t 
times, for some t > 0. Therefore, the length of Q increases by k. □ 

From the above lemmas, we infer that when k units of flow are pushed from v G C to s, 
the length of any path connecting faces adjacent to vertices of C in the dual of ext(P) changes 
by —k, or k. The actual path R v on which we push the flow is not important, the amount 
of change in the length of a path Q is determined only by k, by the parity of the number of 
crossings between Q and R', whose choice was arbitrary and independent of v, and by the cyclic 
order of x, v, and y. 

Moreover, the total amount of change in step 6 in the length of Q is determined only by two 
factors. The absolute value of the change is the total amount of excess that we are returning from 
Ci, it is easy to obtain this value for every pair of dual vertices x, y in total of 0(|C| 2 ) = 0(r) 
time using dynamic programing. The sign of the value is determined by the cyclic order of x, y 
and C\. Note also that all x-to-y paths with a give parity of number of crossing with R' change 
by the same amount, so the shortest paths remain the same. 

We update the lengths of all shortest paths connecting faces adjacent to vertices of C and 
crossing R' an even number of time or odd number of times. Then, for each pair of these faces 
we find smaller of the two distances we have, this way we obtain DDG* (He) for representing 
the sC-flow f[ of value d. Thus, step 6 can be implemented in 0(|C| 2 ) = 0(r) time, since this 
is the number of entries in DDG*(Hc)- 

Determining which side of x and y is C\ is easy, according to the parity of the number of 
crossing between Q and R'. It remains to show how to compute the shortest paths among faces 
adjacent to vertices of C whose number of crossing with R' has a given parity. For this purpose, 
we use a procedure similar to (8j [26] . We build a new graph G R . First, we make an incision in 
ext(P) along R', thus obtaining a graph G R , in which there are two copies of R'\ R\ and i?2- 
Note that the starting and ending vertices of R' are also split. The graph G R is obtained by 
taking two copies of G R and identifying the vertices of R\ in each of the copies of the graph 
with the vertices of i?2 in the other copy of the graph. Each vertex v from the original graph 
has two copies v' and v" in G R . Observe that a path from u' to v" corresponds to a path that 
connects u with v in ext(P) and crosses R' an odd number of times. Similarly, a path from v! to 
v' (or from u" to v") can be mapped to a path crossing R' an even number of times. The same 
is true also if we consider the dual graph of G R , every dual vertex has two copies such that a 
path from one copy to the other corresponds to a dual path in the original graph that crosses 
R' and odd number of times, and a path from a copy back to itself corresponds to a dual path 
in the original graph that crosses R' an even number of times. This means that by using the 
MSSP algorithm of Klein |27j from one of the copies of the dual vertices corresponding to the 
boundary faces along C we can compute the length of the shortest dual paths connecting these 
dual vertices and crossing R' an even or an odd number of times. We uses the length defined by 
the residual capacities with respect to f\. This takes 0((n + |C 2 |)logn) = O(nlogn) time for 
all pairs of dual vertices around C. Note that it is enough to compute these paths only once for 
each piece P, before we apply Algorithm ImplicitMaxFlow. 
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Summing over all steps of the algorithm, we get 0(n log n) preprocessing time per piece, and 
in addition 0(rlog 2 n) running time per sink. Therefore, we obtain the following: 

Lemma 13. Using Algorithm ImplicitMaxFlow we can find an implicit representation of 
a max st-flow, f s t, for every vertex t in a piece different from s, such that f s t is represented 
explicitly in the piece containing t and implicitly for the rest of the graph, in 0(n 2 log n/r + 
nr log 2 n) total time. 

We now show how to use the flow f s t that we found above to store implicitly a min si-cut. 
First note that there are two special cases where this is not required. We already mentioned 
that we explicitly store a min si-cut for every t in the same piece as s using a total of 0{nr) 
space. Another case where we store the si-cut explicitly is when we find that the value fst, for 
some sink t, is the same as the value of /i, the maximum flow from s to the outer boundary of 
the piece containing t. In this case, a min sC-cut is also a min si-cut, there are only 0(n/r) 
different such cuts, and we store all of them explicitly in 0(n 2 /r) space. 

In the remaining case, there is a min st-cut that contains at least one arc of P. We show 
below how to use Lemma [9] to find the cut, which is dual to a separating cycle in the dual planar 
graph. We store the arcs of the cut inside P explicitly, and the arcs of the cut that are not in P 
implicitly, by storing arcs of the dense distance graphs of holes of P that corresponds to these 
parts of the cut. This requires 0(r) space per sink and a total of 0(nr) space. We also keep 
all the MSSP data structures that were used to compute the dense distance dual graph, this 
takes a total of 0(n 2 /r) space. From an MSSP data structure, we can extract a shortest path of 
length t between two vertices in O(^loglogn) time (recall that each arc of a dense distance dual 
graph is defined by a shortest path that is stored in one of the MSSP structures). Note that 
even though the distances in the dual graph that defines the MSSP data structure are changing 
during the execution of the algorithm, as the residual capacity changes, the paths themselves do 
not change. One special case is arc of the dense distance graph that represent the flow in ext(P) 
- in this case the arc corresponds to a shortest path in one of two possible MSSP structures (see 
the details of step 6 of Algorithm ImplicitMaxFlow), we need to store both structures and 
for each arc remember from which of the two it was originated. 

It remains to show how to find the arcs that define the cut. This step is also done during the 
preprocessing of the data structure. Consider the representation of f s t, that is an explicit flow 
assignment for arcs of P and dense distance dual graphs for the holes of P, including ext(P), 
were the distances are defined by the residual capacity with respect to f S f. First, remove all arcs 
of positive length, both in P and in the dense distance dual graphs, since they cannot be in a 
saturated cut. By Lemma [9] we need to find a cycle such that at least one of the arcs of the 
cycle carries flow. For arcs of P, we have the explicit representation of the flow, so it trivial to 
determine if an arc carries a flow. Similarly, for each arc (x, y) in a dense distance dual graph 
of a hole of P, we need to check whether it corresponds to a path that is crossed by some flow. 
Since every arc of the minimum cut has non-negative flow assigned to it, it is enough to compare 
the length of (x, y) defined by the residual capacities with respect to f s t with the length of this 
arc for defined by the original capacity function. 

Thus, the problem can now be stated the following way: given a directed graph with a set 
of distinguished arcs (those who carry flow), find a cycle with at least one distinguished arc. 
We solve this problem in linear time, by first finding strongly connected components and then 
finding a strongly connected component with a distinguished arc (x, y). It follows that there has 
to be a yx path in this strongly connected component. 

The cycle we find consists both of dual arcs of P and arcs from dense distance dual graphs. 
We find the cycle in 0(r) time per sink, and store it using 0{r) space. As we mentioned above, 
on query time we recover the paths inside the holes of P that correspond to these arcs from the 
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MSSP structures in total of 0(|M| log log n) time, where M is the shortest si-separating cycle. 

Recall that in the preliminaries we added some zero-capacity arcs to the graph to make it 
satisfy some structural properties (bounded face size, decomposition with connected pieces). 
Adding these arcs does not change the value of a min si-cut, but might change the size of 
the cut-set. We remove these artificial arcs, and replace all the MSSP data structures with 
equivalent data structures (over the same set of vertices, inside the same subgraph) that store 
shortest paths in the original dual graph. First note that this does not change the distances in 
the MSSP structures (since instead of a dual arc of length we merge the two dual vertices it 
is adjacent to into one). Second, note that a set of dual vertices that are on a single dual face, 
remains on a single dual face also when we remove the artificial arcs (we might need to add 
some infinite length artificial arcs to the dual graph, if the boundary of that dual face becomes 
disconnected). Therefore, by storing the MSSP data structure over the original dual graph we 
can return the cut-sets in time almost proportional to their sizes in the original graph (up to 
the log log n factor we get from using the MSSP data structure). 

This concludes the description of the algorithm, and proves the correctness of the following 
theorem. 

Theorem 2. There exists a data structure that, given a planar n- vertex digraph G = (V,E) with 
a fixed source s in G and a parameter r £ [l,re], can report the min st-cut sets for any sink t. 
The data structure requires 0(n 2 log n/r + rerlog 2 n) time for preprocessing and the min st-cut 
set M is reported in 0(\M\ log log n) time. The data structure requires 0(n 2 /r + nr) space. 

To minimize the required storage space we set r = n 1 / 2 and get 0(n 3 / 2 log 2 n) preprocessing 
time and 0(re 3 / 2 ) space. To minimize the preprocessing time we set r = (n/ log re) 1 / 2 and get 
0(n 3 / 2 log 3 / 2 re) preprocessing time and 0(n 3 / 2 log 1 / 2 n) space. 

If we want to improve the query time to 0(|M|), we can avoid using the MSSP data structure 
for the queries, and store the shortest paths trees it represents explicitly. The MSSP structures 
are defined over total of 0(n/y/r) boundary vertices, each is a member of a constant number of 
structures. We find the explicit shortest path trees from each boundary vertex in linear time per 
tree using the shortest path algorithm of Henzinger et al. [20] . Therefore, the required additional 
time is 0(re 2 j\fr) and the additional storage space is 0(n 2 /y / r) as well. 

Theorem 3. There exists a data structure that, given a planar n-vertex digraph G = (V, E) 
with a fixed source s in G and a parameter r 6 [1, re], can report the min st-cut sets for any sink 
t. The data structure requires 0(n} j ^Jr + nr log 2 re) time for preprocessing and the min st-cut 
set M is reported in 0(\M\) time. The data structure requires 0(n 2 / ' yjr + nr) space. 

To minimize the required storage space we set r = n 2 / 3 and get 0(n 5 / 3 log 2 n) preprocessing 
time with 0(re 5 / 3 ) space. To minimize the preprocessing time we set r = (re/ log 2 re) 2 / 3 and get 
0(n 5 / 3 log 2 / 3 n) preprocessing time with 0(n 5 / 3 log 2 / 3 re) space. 

6 Concluding Remarks 

We gave an 0{n log 3 n) time algorithm for the problem of finding max si- flow values for a fixed 
source s and all sinks t G V \ {s} in an n-vertex planar digraph G = (V,E). The previous best 
known solution was to perform re — 1 executions of a single-source single-sink max flow algorithm, 
which gave an 0(n 2 log n) time bound with the algorithm of Borradaile and Klein [3]. 

An immediate corollary of our result is a near-quadratic time algorithm for finding max 
si- flow values for all source/sink pairs (s, t). We showed that the number of distinct max si-flow 
values is quadratic in the worst case. Hence, our algorithm is optimal up to logarithmic factors. 
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Prior to this result, computing all max- flow values in a directed planar graph required up 
to 0(n 2 ) max-flow computations. We conjecture that a similar improvement is possible for 
general directed graphs. Another interesting direction to pursue is to prove or disprove the 
existence of a data structure similar to the oracle in [6], which can answer queries for max 
si-flow values in constant time after o(n 2 ) preprocessing. A related question is whether it is 
possible to improve the time needed to find max si-flow values for k given input pairs (s,i)? 
Also, is it possible to remove the log-factors, i.e., compute all max-flow values in optimal 0(n 2 ) 
time? Finally, our algorithm computes all cut-sets C in a planar digraph in 0(min{ra 5//2 log 3 ' 2 n+ 
SceC 1^1 1°§ 1°& n ' n§//3 l°g 2 ^ 3 n + SceC l^ll) time. Is it possible to improve this running time? 
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